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We show the possibihty of producing matter-wave switching devices by using Manakov interactions 
between matter wave solitons in two-species Bose-Einstein Condensates (BEC). Our results establish 
the experimental parameters for three interaction regimes in two-species BECs: symmetric and 
asymmetric splitting, down-switching and up-switching. We have studied the dependence upon the 
initial conditions and the kind of interaction between the two components of the BECs. 
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Introduction.- There has been remarkable experimen- 
tal and theoretical progress in recent years regarding the 
interaction phenomena of coherent matter waves within 
Bose-Einstein condensates (BEC) of ultra cold atomic 
gases [H . The interaction between the constitutive bosons 
inside the condensate is defined in terms of the ground 
state scattering length a. When a > the interac- 
tion between the particles in the condensate is repulsive, 
whereas for a < 0, the interaction is attractive. Ex- 
perimental preparation of negative scattering length was 
made possible by using Feshback resonances to contin- 
uously detune the value of a from positive to negative 
values 0. This provides new interest for analyzing sys- 
tems of interacting condensates with an attractive cou- 
pling force, provided the number of particles is limited to 
avoid collapse . 

In spite of this serious difficulty, negative scattering 
length condensates have some peculiarities which make 
them interesting. For instance, if the trap is removed in 
one direction and shrunk in the transverse plane, attrac- 
tive interaction gives rise (for a given number of particles) 
to a self-confined stationary stated . In this case, the 
cloud forms a soliton and can be controlled as a particle 
by acting on it with external fields f5^ l6!|. 

On the other hand, results with two-species BECs0 al- 
low to perform experiments with two-component (vector) 
matter wave solitons. In this configuration, the interac- 
tion between the two species plays a significant role in de- 
termining the dynamics of the clouds |8j. In this work, we 
extend the concept of Bose-Einstein soliton0,|^|^ to the 
case of two-species condensates. We show the possibility 
of producing Manakov solitons^ IT3 | with matter waves 
and to use them in the design of matter-wave switching 
devices. To this aim, we have studied numerically the in- 
teraction dynamics of two-species condensates within the 
confining cylindrical trap for different scattering length 
regimes and initial conditions. 

The theoretical foundations of our work are presented 
in the first part of this paper, where we describe the 
model and the experimental configuration which has been 
assumed. We used a coupled set of nonlinear Schrodinger 



equations, properly describing Manakov solitons [TH. for 
modeling the time evolution of the initial wave functions 
before and after the interaction. Next, we demonstrate 
that some physical insight into the dynamics can be ob- 
tained from a set of equations derived directly from a 
variational analysistu], at least for the case of station- 
ary states, while more complicated scattering phenom- 
ena, must be calculated from direct numerical integration 
of the field equations. We show that interesting switch- 
ing devices can be implemented with these interactions. 
Finally, we present our conclusions. 

Theoretical Foundations.- The system we have stud- 
ied consists of two BEC matter waves of different species 
which interact via 2-body elastic effects. We have as- 
sumed an asymmetric cigar trapping potential in which, 
once individual BEC solitons are formed, they can be 
given an initial momentum along x with an external force 
to make them collide. Assuming that the dynamics of 
the condensate is frozen in the transversal plane y, z and 
the trap is switched off along the x direction, the behav- 
ior of each condensate can be described through a one- 
dimensional Gross-Pitaevskii equation (GPE), with the 
addition of a coupling term that takes into account the 
interaction potential between the two condensate wave 
functions: 
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+ (ofefelV'fcP + ajfc|?/'jp)i/'fe = 0, (I) 



where k = 1,2, j — 2,1. We have defined the dimen- 
sionless variables r = {h/mL'^)t and 77 = x/L, together 
with the normalized wave function (order parameter) 
'iJji{ri,T), which gives the number of particles per unit 
length. Thus, the normalization for tpk is Ni = J \ipi\'^dri, 
being L « lfj,m the radial dimension of the trap, a the 
ground state scattering length (typically a « Sag, with 
ao the Bohr radius) and m the mass of the atoms. 

The analytic solution to the set of equations iQJ rep- 
resents a difficult problem. For some particular cases, 
explicit formulas were given by Manakov f^. However, 
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we can obtain approximate results by assuming an ini- 
tial Gaussian wave function and performing a variational 
analysis This provides physical insight into the prop- 
agation and elastic collisions of fundamental parameters. 
Although this procedure breaks down for inelastic scat- 
tering when the beams split off, we will use these analytic 
results as a guide when we perform a full numerical in- 
tegration of Eqs.QJ. To obtain the differential equations 
for the motion of the centroid and width of the initial 
Gaussian trial wave functions, we must minimize the La- 
grange density, which produces the GP equation, over a 
set of Gaussian trial wave functions such that: 



ipk{'n,T) = Ak - exp 
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, (2) 



where k — 1,2. Inserting this trial function into the La- 
grangian density and taking the variation with respect to 
?7, we obtain a set of differential equations with the above 
parameters, all of them r dependent: Ak (complex am- 
plitude), (T/c (half width of the cloud), rjk (position of 
the centroid), ak (velocity) and Pk (inverse square root 
of the beam curvature radius). The equations obtained 
describe the motion of the centroid and the oscillations 
of the soliton widths. Thus, taking the variation with 
respect to Ak and A^,, and operating with them, we ob- 
tain the conservation of the number of particles: Nk = 0. 
Taking the variation with respect to the parameter ak 
and defining rjjk — rjj — rjk (distance between centroids), 
we obtain the evolution of the separation of the beam 
centers and widths: 
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Eq.||3) can be integrated to obtain the potential 11^ 
ruling the interaction between the two species: 
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with fc = 1,2; j = 2,1. Equations ^ and |Q, although 
not exact, are valuable tools for a further detailed nu- 
merical exploration. In first place, we must notice that 
depending on the sign of the coupling constant, the in- 
teraction potential will correspond to a barrier {ujk < 0) 
or to a well {ajk > 0). Thus, we can predict a mini- 
mum separation distance rjcr as the half width estimate of 



Hk , which corresponds to a value rjcr ~ y'^j ~^ '^k ■ This 
means that two clouds separated more than their width 
will not interact. Our numerical simulations corroborate 
all these predictions from the variational analysis. 

On the other hand, when the two clouds are in the 
same initial position (i.e., 77^(0) — 77fc(0)) with null ve- 
locities, the condition ai = a2 — 0, predict the existence 
of stationary states corresponding to a zero spreading of 
the wave functions. These vector solitons are formed due 
to self-trapping from ajj and akk terms in Eq. |^ to- 
gether with cross-interaction (a^fc). We have performed 
a numerical calculation by means of a relaxation method 
of these stationary states for different values of the scat- 
tering lengths. The results are in agreement with the 
variational calculations. As we will comment below, the 
computer simulations reveal that tunning the value and 
sign of the ajk can yield to a dramatic change in the 
stability of the stationary states. From a physical point 
of view, it is evident that a repulsive cross-interaction 
(ujk < 0) opposed to the self-trapping effect could yield 
to an instabilization of the vector BEG soliton. This situ- 
ation is also predicted by the variational calculation when 
(Tj = in Eq. Q . In this case, the stationary states will 
only exist if: 
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For the particular case of two condensates with equal 
number of particles and widths, the above condition © 
takes the form: \ajk/ajj \ ^ 1. The numerical simulation 
sets this value to \ajk/o-jj\ ^ 0.95, which is in very good 
agreement with the theoretical prediction. 

Numerical simulations. - For the integration of Eqs.(^l, 
we utilized a Crank-Nicholson finite difference discretiza- 
tion with 500 points grid. Our simulations consisted, in 
first place, of comparing different velocities and number 
of particles of the incident BEG soliton wave functions 
described, and observing the evolution of the clouds after 
the interaction, for several scattering lengths. From an 
exhaustive numerical exploration, we realized that there 
exist well defined scattering regimes for which the emerg- 
ing wave functions can be radically different. We also 
explored the effect of tunning the value of the scatter- 
ing length for several collision processes. We have found 
that a change of the sign and the values of the ajj and 
Ojk terms has a deep influence on the dynamics of the 
clouds. 

Thus, we will show in first place the results obtained 
for the case of attractive scattering lengths coeflticients 
(i.e.: ajj and Ojk > 0). Fig. ^ illustrates different Y- 
switching processes based upon the characteristics of the 
incoming wave function. In (a) and (b) we show respec- 
tively the square modulus of ui and U2. The simula- 
tions correspond to the parameters: Ni = N2 = 4000, 
width of the clouds wi = W2 — lO/im, initial separa- 
tion of the centroids xq = 15/im, and initial velocities 
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FIG. 1: [Color on line] Different Y-switching configura- 
tions based upon the characteristics of the incoming wave 
functions. The simulations correspond to the parameters: 
Ni = N2 — 4000, width of the clouds lO/im, initial sep- 
aration of the centroids xo = 15/im, and initial velocities 
Vii — —Vi2 — 0.05/im/s (left), and vn = —Vi2 = 0.5/im/s 
(right). Scattering length parameters are: an — 022 ~ 1.0 
and ai2 = 021 = 2.0. 



Vii — ~Vi2 — 0.05/im/s. The scattcrmg length adimen- 
sional coefficients are an = a22 = 1.0 and ai2 = 021 = 
2.0. The two incoming clouds split symmetrically into 
two mutually trapped vector BEC solitons with the same 
number of particles. In c) the maximum of and 
|w2p are plotted. Figs, ^d and ^e correspond to the 
same parameters but increasing the initial velocities one 
order of magnitude: vn = —Vi2 — 0.5/im/s. In this case 
the emerging distributions are not symmetric, as can be 
clearly appreciated in Fig. ^f where the evolution of 
both peaks of is shown. In both cases the resulting 
matter waves are stable vector BEC solitons with larger 
velocities than the input beams. The wave functions take 
the form aiUi -|- a2M2, being |aip -I- \oi2\^ — 1. We have 
found numerically that the different splitting regimes are 
limited by the following values: < \vi\ < O.l/im/s for 
the symmetric case and O.l/zm/s < \vi\ < 0.9/im/s, for 
the splitting with different number of particles. For ini- 
tial velocities \vi\ > 0.9/im/s the two wave functions are 
mutually transparent and arise with any change after the 
collision. 

The collision of two vector BEC solitons constructed 
by superposition of Gaussian-shaped clouds is shown in 
Figure |21 In this case it is obtained an up-switching 
of the peak of both wave functions which emerge fused 
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FIG. 2: [Color on line] Up-switching of two condensate pairs 
with different number of particles. A''i; — Nir = 8000 (a) 
and N21 = N2r ~ 3000 (b); widths 5/im; initial velocities 
vi = —V2 — O.Sfim/s and initial separation of the centroids 
Xo = 15/im. Scattering length coefficients: an = a22 = 1.0 
and ai2 — 021 = 2.0. (c) and (d) show the evolution of the 
peaks of (a) and (b), respectively. Insets show input and 
output wavefuntions for the respective cases. 



after the collision. The simulation corresponds to the 
parameters: number of particles: A^i; = iVi^ — 8000 
(a) and N21 = N2r = 3000, where subindices /, r des- 
ignate the soliton coming form left or right side (b); 
cloud widths 5/im; initial velocities vi — — U2 = 0.3fim/s 
and initial separation of the centroids a;o — 15/im. The 
scattering length coefficients are: an — 022 = 1.0 and 
ai2 = 021 = 2.0. We must stress that the fusion of 
two condensates only takes place for clouds with different 
number of particles. 

We will now consider the case of repulsive cross- 
interaction {ajk < 0). This situation is very interest- 
ing, as it corresponds to condensates of different types of 
atoms (for instance ^^Rb and '^'^K) as in the experiments 
from ref.Q. As it can be appreciated in Fig. |2| the results 
of the collisions are very different comparing with Fig. ^ 
which corresponds to attractive cross-interaction. Figs. 13 
(a) and (b) show the down-switching of one of the conden- 
sates which is almost annihilated by the other, which is 
repelled. The plot (c) shows the evolution of the peaks of 
both distributions. The parameters used for the calcula- 
tion are typical of experiments: Ni = N2 = 4000, widths 
of the clouds b.Ofim, initial separation of the centroids 
xq = 30/im, and input velocities vi — —V2 — Q.S^m/ s. 
The values of adimensional scattering length parameters 
an — 0.6, 022 = 1.0, and ai2 — 021 — —2.3 were taken 
form ref. Q . In figures|31-d and|31-f we show the same sim- 
ulation but taking an — a22 = 1-0, and ai2 — 021 = — 1 
In this case the effect of the collision is dramatic and both 
clouds almost completely spread. 

Finally we have studied the evolution of two initially 
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FIG. 3: [Color on line] Down-switching of one condensate 
(left) and two condensates (right), in the case of repulsive 
cross-interaction. The simulations correspond to the parame- 
ters: number of particles A''i — N2 ~ 4000, width of the clouds 
10.0/im, initial separation of the centroids xq — SO.Ofim, and 
initial velocities vi — V2 = 0.8fim/s. The scattering length 
parameters are: an = 0.6, 022 = 1.0, and ai2 — 021 ~ —2.3 
(left); ai2 — ^21 = —1 and an = a22 = 1.0 (right). Peak 
power plots of the condensates before and after the collision 
are shown in (c) and (f) for respective cases. 

superposed condensates in the case of repulsive cross- 
interaction, with Oil = 022 — 1.0 and ai2 = 021 = —0.8. 
The widths of the clouds are lO^m and Ni ^ N2 = 15500 
have been chosen to satisfy the condition of stationary 
state from Eq. We have observed that for clouds 



with almost exact number of particles, the existence of 
the stationary state is not affected by the negative sign 
of ajk. However, if the number of particles is slightly al- 
tered (taking, for instance Ni — 15430 and N2 — 15750), 
the condensate with less population splits off in two soli- 
tons and the other condensate remains oscillating, as 
shown in Fig^l 

Conclusions.- We have studied, through an approxi- 
mate theoretical derivation and through numerical sim- 
ulations, several consequences of the interaction of com- 
posite Bose Einstein mutually interacting solitons. These 
systems give rise to a mutual trapping phenomena due 
to the soliton-soliton interaction forces. As a result of 
these interactions, we have predicted new phenomena 
not previously detected for the formation of stable vector 
matter-wave solitons, after a violent division, and have 
shown that there exist a critical initial velocity for these 
(a) , ■ (b) 
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FIG. 4: [Color on line] Propagation of two initially super- 
posed condensates in the case of repulsive cross-interaction, 
with a slightly different number of particles A'^i = 15430 
and N2 = 15750. The widths of the clouds are lO/im and 
the scattering length parameters are an = a22 ~ 1.0 and 
ai2 = "21 = —0.8. 



phenomena to occur. Depending on the values and signs 
of the scattering length coefficients, we have found sev- 
eral different switching regimes for matter wave solitons. 
Finally, we have calculated the effect of a repulsive cross- 
interaction in the stability of vector Bose-Einstein soli- 
tons. 
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